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Abstract. This paper studies the reaction of low-mass 
stars to anisotropic irradiation and its implications for the 
long-term evolution of compact binaries (cataclysmic vari- 
ables and low-mass X-ray binaries). 

First, we show by means of simple homology consid- 
erations that if the energy outflow through the surface 
layers of a low-mass main sequence star is blocked over a 
fraction s e g < 1 of its surface (e.g. as a consequence of 
anisotropic irradiation) it will inflate only modestly, by a 
factor ~ (1 — Sag) . The maximum contribution to mass 
transfer of the thermal relaxation of the donor star is s e g 
times what one obtains for isotropic (s e g = 1) irradiation. 
The duration of this irradiation-enhanced mass transfer 
is of the order of 0.1 1 ln(l — s e g)| times the thermal time 
scale of the convective envelope. Numerical computations 
involving full ID stellar models confirm these results. 

Second, we present a simple analytic one-zone model 
for computing the blocking effect by irradiation which 
gives results in acceptable quantitative agreement with 
detailed numerical computations. 

Third, we show in a detailed stability analysis that 
if mass transfer is not strongly enhanced by consequential 
angular momentum losses, cataclysmic variables are stable 
against irradiation-induced runaway mass transfer if the 
mass of the main sequence donor is M ^ O.7M0. If M <; 
O.7M systems may be unstable, subject to the efficiency 
of irradiation. Low-mass X-ray binaries, despite providing 
much higher irradiating fluxes, are even less susceptible to 
this instability. 

If a binary is unstable, mass transfer must evolve 
through a limit cycle in which phases of irradiation- 
induced high mass transfer alternate with phases of small 
(or no) mass transfer. At the peak rate mass transfer pro- 
ceeds on s ff times the thermal time scale rate of the con- 
vective envelope. A necessary condition for the cycles to be 
maintained is that this time scale has to be much shorter 
(^ 0.05) than the time scale on which mass transfer is 
driven. 



Send offprint requests to: hsr@mpa-garching.mpg.de 



Key words: accretion, accretion disks — binaries: close 
- Stars: evolution — Stars: interiors — novae, cata- 
clysmic variables X-rays: stars 



1. Introduction 

It is only a few years since the realization that irradia- 
tion of the donor star of a semi-detached binary by ac- 
cretion luminosity generated in the vicinity of a com- 
pact accretor can have far-reaching consequences for the 
long-term evolution of such binary systems. Podsiadlowski 
(1991), treating irradiation in spherical symmetry, first 
addressed this problem in the context of the evolution of 
low-mass X-ray binaries (LMXBs). Later, Gontikakis and 
Hameury (1993), Hameury et al. (1993), and Ritter (1994) 
pointed out that spherically symmetric irradiation is prob- 
ably not an adequate approximation for the case in point. 
Preliminary calculations by Ritter (1994) showed that 
the reaction of low-mass stars to anisotropic irradiation 
is qualitatively (and quantitatively) different from that 
to isotropic irradiation. Subsequently, King et al. (1995, 
1996, 1997), hereafter KFKR95, KFKR96 and KFKR97, 
Ritter, Zhang and Kolb (1995, 1996), Ritter, Zhang and 
Hameury (1996), Hameury and Ritter (1997), hereafter 
HR97, and McCormick and Frank (1998), hereafter MF98, 
have dealt with the case of anisotropic irradiation in var- 
ious ways and detail. In most of the above-quoted papers 
reference was made to a paper by the present authors in 
which a detailed and systematic treatment of the basic 
properties of anisotropic irradiation of low-mass stars was 
to be found. The preparation of this material has been 
much delayed. The main purpose of the present paper is 
to close this gap in the recent literature, and to provide 
the base for further research. 

For this purpose we shall discuss in section 2 in some 
detail the basic concepts of computing the mass transfer 
rate in a semi-detached binary with and without irradi- 
ation, and present observational evidence and theoretical 
arguments in support of the notion of anisotropic irra- 
diation. In section 3 we shall present a simple analytic 
calculation which shows that the reaction of a low-mass 
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star to anisotropic irradiation is qualitatively and quanti- 
tatively different from what one obtains in the isotropic 
case. We show also that the simple analytic results are 
fully supported by corresponding calculations of full ID 
stellar models (hereafter simply referred to as full stellar 
models) . In section 4 we study the stability of mass trans- 
fer in binaries in which the donor star is exposed to irradi- 
ation generated by the accreting companion. In section 5 
we shall present a more detailed semi-analytic model than 
the one in Sect. 3 for the reaction of a low- mass star to 
anisotropic irradiation in the limit of small fluxes. The 
implications of the instability against irradiation-induced 
runaway mass transfer for the long-term evolution of CVs 
and LMXBs is discussed in section 6. In section 7 we shall 
present and discuss results of numerical computations of 
the secular evolution of CVs subject to the irradiation in- 
stability. A summary of our main results and our main 
conclusions are given in the final section 8. 



2. Computing the mass transfer rate 

In the context of this paper, compact binaries are either 
cataclysmic variables (CVs) or LMXBs, i.e. systems con- 
sisting of a low-mass star (the secondary) with a mass 
M s ^ 1Mq which fills its critical Roche lobe and transfers 
matter to a compact companion (the primary), either a 
white dwarf (in CVs), or a neutron star or a black hole 
(in LMXBs). The secular evolution of such systems is a 
consequence of mass loss from the secondary which, in 
turn, is driven either by the secondary's nuclear evolu- 
tion or by loss of orbital angular momentum and possibly 
other mechanisms such as irradiation on which this paper 
focuses. In the standard picture of the secular evolution 
see (e.g. King 1988; Kolb and Ritter 1992, hereafter KR92; 
Ritter 1996) the nature of the compact star is of no impor- 
tance, i.e. the star is considered as a point mass (of mass 
M c ). The nature of the compact component is, however, 
of importance for accretion phenomena occurring in such 
systems, e.g. dwarf nova and classical nova outbursts in 
CVs and X-ray bursts in LMXBs, and when illumination 
of the secondary by radiation generated through accretion 
is considered, as we shall do in the following. Thus, cal- 
culating the standard secular evolution of such a binary 
reduces to calculating the evolution of a low-mass star un- 
der mass loss, where the mass loss rate derives from the 
boundary conditions imposed by the fact that the star is 
in a binary. In the simplest case one obtains the mass loss 
rate M s by requiring that the radius of the secondary R s 
is always exactly equal to its critical Roche radius Rr. By 
decomposing the temporal change of R s and Rr as (see 
e.g. Ritter 1988, 1996) 
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is the adiabatic mass radius exponent of the secondary, 
and 



Cr 



dlnR.R 
SlnAG 



(5) 



the mass radius exponent of the critical Roche radius, 
where the subscript * indicates that for evaluating this 
quantity one has to specify how mass and angular momen- 
tum are redistributed in the system. J is the orbital an- 
gular momentum, (d\nR s /dt) t h the rate of change of R s 
due to thermal relaxation and (dlaR s /dt) nuc the one due 
to nuclear evolution. The virtue of (||) is that it shows im- 
mediately how secular evolution works: If the binary is dy- 
namically stable against mass transfer; i.e. if Cs — Cr > 0, 
then mass transfer must be driven by some mechanism. 
This can either be the secondary's expansion due to nu- 
clear evolution or the shrinkage of the orbit due to loss of 
orbital angular momentum. In the standard evolution of 
low-mass binaries, mass transfer is usually not driven by 
thermal relaxation. However, as we shall see below, this is 
not necessarily true if irradiation of the secondary is taken 
into account. 

Attempts to bring the observed population of millisec- 
ond pulsars in line with the death rate of their presumed 
progenitors, i.e. the LMXBs, have resulted in the spec- 
ulation that the secular evolution of LMXBs might be 
drastically different from, and their lifetime much shorter 
than that of CVs (e.g. Kulkarni and Narayan 1988). A 
possible reason for this is seen in the fact that the donor 
in a LMXB is exposed to a high flux of hard X-ray radia- 
tion emitted by the accreting compact star. In fact, Pod- 
siadlowski (1991) has shown that irradiating a low-mass 
main-sequence star (M ^ 0.8M Q ) spherically symmetri- 
cally with a flux Fi rr J> 10 11 — 10 12 erg cm _2 s _1 results in 
an expansion of the star on a thermal time scale and in 
its gradual transformation into a fully radiative star. It is 
this expansion which can drive mass transfer on a thermal 
time scale and this was thought to shorten the lifetime of 
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a LMXB. On the formal level the effect of irradiation is 
taken into account by replacing ([!]) by 



dlnR s _ M s 

dt ~ ^ s mZ 
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Here, (d\nR s /dt) m \ is the thermal relaxation term aris- 
ing from mass loss alone and (3 In R s /dt)i rl the thermal 
relaxation term caused by irradiation (at constant mass). 
The latter term arises because irradiating the donor star 
means that its surface boundary conditions are changed 
and that it tries to adjust to them on a thermal time scale. 
With (|) instead of (g) we have now for the mass transfer 
rate 
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From (Q) it is clearly seen that if (3 In R s /dt) m > 0, irra- 
diation amplifies mass transfer. However, there are limits 
to what irradiation can do because the time scale and the 
amplitude of the effect are limited. The secondary's expan- 
sion due to irradiation saturates at the latest when it has 
become fully radiative, where it is larger by 5hiR s than a 
star of the same mass in thermal equilibrium without irra- 
diation. As we shall show below the time scale r on which 
the irradiated star initially expands is of the order of or 
less than the thermal time scale of the convective envelope 
t cc . As a result, we can expect a maximum contribution 
from irradiation to M s at the onset of irradiation (i.e. of 
mass transfer) which is of the order 



(-6M s ) in = 



Cs-b 



dt 
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However, with ongoing irradiation not only does 
{d\nR s /dt)- m decrease roughly exponentially on the time 
scale t, but other effects also tend to decrease — M s . With 
the star's transformation from a mainly convective to a 
fully radiative structure, Cs increases from — 1/3 ^ Cs <0 
to a large positive value. Thus Cs — Cr increases, with 
the result that —M s decreases. In addition, even the driv- 
ing angular momentum loss is affected: for both brak- 
ing mechanisms discussed in the literature, i.e. magnetic 
braking (e.g. Verbunt and Zwaan 1981; Mestel and Spruit 
1987) and gravitational radiation (e.g. Kraft et al. 1962), 
(9 In | J\/ d\nR s ) < 0. Thus \J\ decreases if the secondary 
inflates. Worse, magnetic braking which is thought to be 
coupled to the presence of a convective envelope might 
cease altogether once the star has become fully radia- 
tive. In this case the system might no longer be able 
to sustain the mass transfer necessary to keep the sec- 
ondary in its fully radiative state. As a consequence, it 
is therefore possible that such systems evolve through a 



limit cycle in which short phases (of duration r) with 
irradiation-enhanced mass transfer alternate with long 
detached phases during which the oversized secondary 
shrinks as it approaches thermal equilibrium, again on the 
thermal time scale, whereas the system contracts on the 
much longer time scale of angular momentum loss. 

A number of theoretical studies along these lines, un- 
der the assumption of spherically symmetric irradiation, 
have been performed (e.g. by Podsiadlowski 1991; Harpaz 
and Rappaport 1991; Frank et al. 1992; D'Antona and 
Ergma 1993; Hameury et al. 1993; Vilhu et al. 1994), all 
more or less confirming the behaviour outlined above, in- 
cluding cyclic evolution (Hameury et al. 1993; Vilhu et al. 
1994). It should be stressed once more that the validity 
of these studies, in which the effect of irradiation is maxi- 
mized, rests on the validity of the assumption that spheri- 
cally symmetric irradiation is an adequate approximation. 
In fact, Gontikakis and Hameury (1993) and Hameury et 
al. (1993) have examined whether the spherically sym- 
metric approximation is adequate and find that it is not. 
Worse, when taking into account the anisotropy of the 
irradiation they find that the long-term evolution differs 
significantly from the spherically symmetric case. Now, 
anisotropic irradiation is a rather difficult 3-dimensional 
problem involving hydrostatic disequilibrium to some ex- 
tent, and with it circulations which can transport heat 
from the hot to the cool side of the star. Because of this one 
might dismiss simple theoretical arguments such as those 
given by Gontikakis and Hameury (1993) and Hameury 
et al. (1993). There is, however, a much more compelling 
argument supporting the case of anisotropic irradiation. 
This derives from the observations of a number of close but 
detached binary systems in which a low-mass companion, 
exposed to an intense radiation field emerging from its 
hot (degenerate) companion, shows a bright illuminated 
and an essentially undisturbed cool hemisphere. A com- 
pilation of the systems in question is given in Table 1, 
where we list the object's name, its association with a 
known planetary nebula, its orbital period (in days), the 
effective temperature of the irradiating white dwarf, the 
amplitude (in magnitudes) of the "reflection effect", and 
relevant references. These systems all show a pronounced 
"reflection effect" in their light curves which, in turn, is ex- 
plained by the anisotropic temperature distribution on the 
irradiated companion. Among them are 7 binary central 
stars of planetary nebulae. The systems listed in Table 1 
demonstrate that a cool star, exposed to strong irradi- 
ation from a hot companion, can live without problems 
with a hot and a cool hemisphere. Moreover equal effec- 
tive temperature over the whole surface is not established, 
at least not over atime scale ~ 10 4 yr associated with the 
age of central stars of planetary nebulae, despite the fact 
that the irradiated star need not even rotate nearly syn- 
chronously. Obviously, heat transport from the hot side 
to the cool side by means of circulations is sufficiently 
ineffective that the large difference in effective tempera- 
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Table 1. Detached close binaries showing a reflection effect 



Object 


associated 
planetary 
nebula 


Porb(d) 


T eS ,l{W 3 K) 


amplitude of 
reflection 
effect 


Ref. 


BE UMa 


+ 


2.2912 


105 ± 1 


Am pg = 0.9 


1.2 


VW Pyx 


Kl-2 


0.6701 


85 ±6 


Am = 1.5 


3,4 


V664 Cas 


HFG 1 


0.5817 




AB = 1.1 


5,6 


V477 Lyr 


Abell 46 


0.4717 


60 ±10 


AV = 0.6 


7 


UU Sge 


Abell 63 


0.4651 


117.5 ± 12.5 


AV = 0.4 


8,9 


KV Vel 


DS 1 


0.3571 


77 ±3 


AV = 0.56 


10,11 


TW Crv 




0.3276 




AV = 0.85 


12 


AA Dor 




0.2615 


40 ±3 


AV = 0.05 


11,13,14 


MS Peg 




0.1737 


22.2 ± 0.1 


AV = 0.1 


15 


NN Ser 




0.1301 


55 ±8 


AV = 0.4 


16,17 


HW Vir 




0.1167 


33.0 ± 0.8 


AV = 0.2 


11,18,19 


MT Ser 


Abell 41 


0.1132 


80 ±20 


AB = 0.15 


20,21 


NY Vir 




0.1010 


33 ±3 


AV = 0.2 


22 



References: 1: Ferguson et al. (1999); 2: Wood et al. (1995); 3: Kohoutek & Schnur 



(1982); 4: Bond & Grauer (1987); 5: Grauer et al. (1987); 6: Acker & Stenholm 
(1990); 7: Pollacco & Bell (1994); 8: Pollacco & Bell (1993); 9: Bell et al. (1994); 
10: Landolt & Drilling (1986); 11: Hilditch et al. (1996); 12: Chen et al. (1995); 
13: Kilkenny et al. (1979); 14: Kudritzki et al. (1982); 15: Schmidt et al. (1994); 
16: Haefner (1989); 17: Catalan et al. (1994); 18: Wlodarczyk & Olszewski (1994); 
19: Wood & Saffer (1999); 20: Grauer & Bond (1983); 21: Grauer (1985); 
22: Kilkenny et al. (1998). 



tures can be maintained over long times. Thus the case of 
anisotropic irradiation has to be taken seriously and de- 
serves a more detailed study. This is the objective of this 
paper. 

3. Reaction of a low-mass star upon reducing its 
effective surface 

We start our examination of anisotropic irradiation by re- 
calling one of the main results obtained in the studies 
of spherically symmetric irradiation (e.g. Podsiadlowski 
1991; D'Antona and Ergma 1993), namely that the main 
effect of irradiating a low-mass main-sequence star is that 
the star cannot lose energy as effectively (or at all) through 
the irradiated parts of its surface. If irradiation is spher- 
ically symmetric, the star has no choice but to store 
the blocked energy in gravitational and internal energy 
with the well-known result that it swells. In the case of 
anisotropic irradiation, such as one-sided irradiation from 
an accreting companion, the situation is qualitatively dif- 
ferent: in addition to storing the blocked luminosity in 
internal and gravitational energy, the star can also di- 
vert its energy flow to and lose energy from the unir- 
radiated parts of its surface. This is easily possible in 
the adiabatic convection zone. Because in this zone the 
energy flow is almost fully decoupled from the mechani- 
cal and thermal structure (the flow being proportional to 
(V - V Q ) 3/2 with V-V««l rather than to V (where 
V = d In T/d In P is the actual temperature gradient and 
V a = (91nT/<91nP) a the adiabatic temperature gradient, 
and T and P are respectively the temperature and pres- 



sure) the mechanical and thermal structure of the star can 
still be considered to be spherically symmetric despite the 
fact that the energy flow might be highly anisotropic. The 
only thing which changes is the surface boundary condi- 
tion which replaces the Stefan-Boltzmann law in the unir- 
radiated case. 

In the following we shall make a simple model for 
studying the situation described above. In this model we 
assume that over a fraction s c g of the stellar surface the 
energy outflow is totally blocked (because of irradiation, 
or e.g. star spots, see Spruit and Ritter 1983) and that the 
remaining fraction of the surface (1 — s cff ) radiates with an 
effective temperature T e g. The surface luminosity of the 
star L can therefore be written as 

L(s eff ) = AttR 2 {\ - s eS )aT* s , < s cS < 1 , (9) 

where R is the stellar radius. It is clear that in a more 
realistic model s e ff itself must depend on the irradiating 
flux F{ TT . In section 4 we shall discuss results of numeri- 
cal calculations and a model with which we can evaluate 
Seff(-Flrr)- For the moment we note only that s c ff — > as 
Fi IT — ► and that in the limit of high F{„, s e ff approaches 
the surface fraction of the star which sees the irradiation 
source. 

The purpose of this simple modelling is, on the one 
hand, to provide estimates for the magnitude of and the 
time scale associated with the thermal relaxation pro- 
cess enforced by anisotropic irradiation, and, on the other 
hand, to show that the effects of anisotropic irradiation 
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are not only quantitatively but also qualitatively different 
from those obtained in the spherically symmetric case. 

The internal structure of a low-mass star with a deep 
outer convection zone can be described by the simple ana- 
lytical model by Kippenhahn and Weigert (1994) for stars 
on or near the Hayashi line. This model is particularly 
applicable in our case, since we are only interested in the 
differential behaviour. Assuming a power law approxima- 
tion for the frequency independent, i.e. grey photospheric 
opacity « p h on the unirradiated surface of the form 



Kph = const. P a T b 



(10) 



the Eddington approximation yields the photospheric so- 
lution 

logT cff = -^-ilogP ph + ilogA/-|logi?+const.,(ll) 

where P p h = P(r = 2/3) is the photospheric pressure at 
an optical depth r = 2/3. The interior solution can be 
approximated by a polytrope with index n = 3/2 and 
yields 



logT= | logP+^ 



3 1 

- log R — — log M + const. 



(12) 



Taking ( |12| ) at the photospheric point (P = P p \ r , T = 
T e g), i.e. equating ( |TT| ) and (p^), together with (||) and 
( |Io| ) yields the luminosity L lost by the star with radius 
R over the surface area 47ri? 2 (l — s c g): 



L(s eS ) = L (l - s cS ) 



R 

Rq 



(22a+4&+6)/(5a+2h+5) 



(13) 



Here Lq and Rq are respectively the luminosity and the 
radius of the unirradiated star in thermal equilibrium. In 
general, the irradiated star will not be in thermal equilib- 
rium, i.e. L(R, s e ff) does not equal the nuclear luminosity 
L nuc . For the disturbed star, the latter can be estimated 
by using homology relations. If we write the rate of nu- 
clear energy generation e nuc in the form appropriate for 
hydrogen burning via the pp-chain, i.e. 



£nuc = const. qT v , 



(14) 



where g is the density, we obtain (see e.g. Kippenhahn and 
Weigert 1994) 



Lnnc(R) — Lq ( — 
Kq 



(15) 



From (O) and (fT5J) we obtain the gravitational luminosity 



(16) 



L g (R, s cff ) — L(R, s cff ) - L nuc (R) . 



L g (R,s c ff) — defines the thermal equilibrium values of 
the irradiated star. These can be expressed in terms of 



the thermal equilibrium values of the unirradiated star as 
follows: 



Re(s e ff) = i? (l - S e ff) r 



with 



(17) 



5a + 2b + 5 



(^ + 3)(5a + 26 + 5) + 22a + 46 + 6 ' 



with 



and 



with 



L e (s e g) = Lq(1 - S e s) 



Tes(Sefi) = r e ff,o(l — Scff)* 



(18) 



(19) 



t = 



l-2r-\ 



Since the dominant opacity source in the photosphere of 
cool stars is H~ bound-free absorption, appropriate values 
for a and b in ( |l(]| ) are a sa 0.5 and b w 4—5. Nuclear energy 
generation in low-mass stars occurs mainly via the pp-I 
chain and the appropriate value of v in ( |l4| ) is v w 3 — 5. 
As a result, we find that 

r w -0.1 



and 



t 



0.75 



-0.003 



This means that the effective temperature on the unir- 
radiated part of the surface hardly changes, reflecting a 
well-known property of stars on or near the Hayashi line. 
Furthermore, since r w —0.1, the response of the stel- 
lar radius to anisotropic irradiation (s e g < 1) is much 
weaker than if isotropic irradiation is assumed. Specifi- 
cally, if s e ff ^ 0.5, as is the case for one-sided irradiation 
by an accreting companion, the equilibrium radius with 
irradiation is larger than the one without irradiation by 
at most ~ 7%. On the other hand, the total luminosity L e 
is significantly reduced: because of the rather large value 
ofi/ + 3~6 — 8 already a slight expansion of the star leads 
to a marked reduction of its nuclear luminosity. 

Clearly, our model is not applicable if s e ff — > 1 ■ This is 
because on the formal level L — > if s c g — > 1 (see Eq. |9|) 
and the values of R and T c s in thermal equilibrium with 
r < and t < (Eqs. [l7] and |l^) diverge. There is also 
a physical reason why this model does not apply in this 
case: s c g — 1 corresponds to strong, spherically symmetric 
irradiation where the star in thermal equilibrium is fully 
radiative (e.g. Podsiadlowski 1991), whereas our model 
applies only to the extent that the star in question retains 
a deep outer convective envelope, even when irradiated. 

We can now estimate the maximum contribution to 
mass transfer arising from thermal relaxation due to irra- 
diation and the associated time scale. Thermal relaxation 
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Fig. 1. The function T(Q,ni) defined in Eq. @ (full 
line) and J-(Q, ni)M ce /M (dashed line) for age-zero main 
sequence models in the bipolytrope approximation as a 
function of mass. 

is maximal right at the onset of irradiation. At that time 
the gravitational luminosity of the star is L g (Ro, s e ff) = 
— s c sLq. Using now the bipolytrope model in the formu- 
lation of KR92, we can write for the thermal relaxation 
term due to irradiation 



d In R 
M 



R=Ro 



GM 2 



Seff nr^i \ 

Tee 



(20) 



where tkh = GM 2 / RqLq is the Kelvin-Helmholtz time of 
the unirradiated star in thermal equilibrium, t cc = tkh jT 
the thermal time scale of the convective envelope, and the 
quantity T{Q,n\) defines a dimcnsionlcss number which 
depends only on the relative size Q = R CO i C /Rq of and on 
the polytropic index n\ in the radiative core. An explicit 
expression for T{Q,ni) is given in KFKR96). In particu- 
lar, for a single polytrope n = 3/2, i.e. a fully convective 
star, where Q = 0, one has T = 7/3. Furthermore, we 
show J-{Q, nx) as a function of mass for zero age main se- 
quence stars as a full line in Fig. 1. We note that T scales 
roughly as the inverse of the relative mass M ce /M of the 
convective envelope: the dashed line in Fig. 1 shows that 
J- ■ M ce /M ss const, within better than a factor of two. 
Thus for the purpose of an estimate we can rewrite ( pp[ ) 
as 



'dlniT 



R=Ro 



7 S e ff (M ce 
3 T KH \ M 



(21) 



Since the radius in thermal equilibrium with irradiation is 
larger by 



SlnR = In 



-Re(Scff) 

Ro 



rln(l - s e ff) 



(22) 



the time scale for thermal relaxation becomes 



6\nR 



dlnR 
dt 



ln(l - s cff ) 

r T r , 



Scff 



ln(l - s off ) 7 M cc 
< r r T KH ~YT 

Seff 3 M 



(23) 



As can be seen from (21) the maximal rate of expansion 
is proportional to s c ff- Interestingly, the time over which 
the new thermal equilibrium is established is much shorter 
than if s c a = 1, unless 1 — s ff is a small number and our 
model does not apply anyway. 

Either (|l]) or ( p3| ) are to be inserted in (||) to get 
an estimate for the contribution of irradiation to mass 
transfer. Using ( pl| ) in (@) yields the peak contribution 
and (|2^) in (||) a time average. 

In order to check the validity of our simple analytical 
model we have also performed numerical computations of 
full stellar models using the modified Stefan-Boltzmann 
law as one of the outer boundary conditions and with 
s c ff as a free parameter. For our computations we have 
used a modified version of Mazzitelli's (1989) stellar evo- 
lution code which is described in more detail in KR92, 
and have assumed a standard Pop. I chemical composi- 
tion with (in the usual notation) X = 0.70 and Z — 0.02. 

One of the basic predictions of our simple model is 
that (see Eq. |j]) log(i? e (s e ff)/i?o) scales linearly with 
log(l — s e ff) and that the slope \r\ is small. This is nicely 
confirmed by the behaviour of full stellar models shown 
in Fig. 2. As can be seen, the prediction is valid, at least 
qualitatively, over more than two orders of magnitude of 
(1 — s c ff). Furthermore the slope in Fig. 2 is indeed small, 
confirming that \r\ is a small number. That the slope is 
different for stars of different mass is due to the fact that 
the effective values of the parameters a, b and v change 
with stellar mass. 

In Fig. 3 we show in a mass radius diagram the ther- 
mal equilibrium radius R e as a function of mass for three 
different values of s c ff, i.e. for the standard main sequence 
(s e ff = 0), and for s e ff = 0.5 and s e s = 0.9. As can be seen, 
the mass radius relations for s e g = 0.5 and s e g = 0.9 are 
shifted by a small amount to larger radii and run roughly 
"parallel" to the standard main sequence (s e ff = 0), as is 
predicted by our simple model. 

Finally, in Fig. 4 we show as an example the thermal 
relaxation with time of an O.4M star with s ff = 0.5. If 
time is measured in units of {dt/ dhiR)\R = R , as is done 
in Fig. 4, it is seen that the relaxation process is character- 
ized by this time scale and that new thermal equilibrium 
is reached after only 0.2-0.3 of these time units. Again 
this confirms our analytical result (cf. Eq. ^3|), according 
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? -2.5 -2 -1.5 -1 -0.5 




-1 -0.5 

lo § (1 -s.fi) 

Fig. 2. Radius i? e (s c ff) of homogeneous stellar models in 
thermal equilibrium as a function of the fraction s c s of 
the stellar surface over which energy loss is assumed to be 
blocked. The full lines in both panels show the prediction 
of the homology model (Eq. [Tt]) with the parameters a = 
0.5, b — 4, and v as indicated. The dashed line in the upper 
panel is for full stellar models with a mass of O.4M0, the 
dashed and dotted lines in the lower panel for full stellar 
models with a mass of respectively O.2M and O.8M . 



to which the relaxation process lasts a few rln(l — s c s) 
in these units. It is also seen that the effective tempera- 
ture on the unirradiated part of the star rises only very 
little, as predicted, but that the relative mass of the con- 
vective envelope is reduced significantly from ~ 0.6 at the 
beginning to ~ 0.33 in the new thermal equilibrium. 

4. Stability against irradiation-induced mass 
transfer 

Let us now examine the situation in which a low-mass 
star transfers mass to a compact companion (of mass M c 
and radius R c ) and, in turn, is irradiated (directly or in- 
directly) by the accretion light source. Because irradia- 
tion can enhance mass transfer and more irradiation can 
give rise to even higher mass transfer, we must examine 
under which conditions such a situation is stable against 
irradiation-induced runaway mass transfer. 

4-1. Arbitrary irradiation geometry 

In this subsection we wish to keep the discussion as gen- 
eral as possible. Therefore, we do not specify the irradia- 




-1 -0.5 

log M/M 

Fig. 3. Mass radius diagram of full stellar models in ther- 
mal equilibrium forr three different values of the fraction 
s c tf of the stellar surface over which energy loss is as- 
sumed to be blocked. Full line: s c s — (normal main 
sequence); short dashed line: s g = 0.5, and long dashed 
line: s e g = 0.9. 

tion geometry. Specific models which do just that will be 
presented in the next subsection (4.2). 

The component normal to the stellar surface of the 
irradiating flux generated by accretion can be written as 

47r K c or 

Here — M s is the mass transfer rate, a the orbital sep- 
aration, ip) a dimensionless function of the position 
on the secondary's surface (characterized by polar coor- 
dinates cp)) which describes the irradiation geometry, 
and rj < 1 a dimensionless efficiency factor, rj = 1 only 
if the accretion luminosity is radiated isotropically from 
the compact star and if all the energy which is radiated 
into the solid angle subtended by the donor as seen from 
the accretor is absorbed below the photosphere. In a real 
situation r\ will be considerably less than unity for several 
reasons. The most important of these are: 1) the accre- 
tion luminosity will, in general, not be emitted isotropi- 
cally, e.g. if accretion occurs via a disk which radiates pre- 
dominantly perpendicular to the orbital plane and casts 
a shadow onto the donor, or if the accretor is strongly 
magnetized and accretes mainly near the magnetic poles, 
as e.g. in AM Her systems. 2) energy emitted in certain 
spectral ranges, as e.g. EUV radiation and soft X-rays, 
is unlikely to reach the photosphere of the donor (in the 
case of EUV and soft X-ray radiation because of the high 
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Fig. 4. Thermal relaxation of a 0.4M Q main sequence star 
after blocking the energy loss over a fraction s e a = 0.5 of 
its surface at time t = to. Top frame: radius R, second 
frame: luminosity L, third frame: effective temperature 
T e ff of the radiating part, and bottom frame: the relative 
mass M ce /M of the convective envelope as a function of 
time. Time is measured in units of the time scale on which 
the radius grows at t = to- i?o, Lq, and T e ff,o are respec- 
tively the values of R, L and T cr f immediately before the 
onset of the blocking of energy outflow. 



column density of neutral hydrogen). 3) part of the inci- 
dent flux will be scattered away before penetrating into 
the photosphere. 4) not all of the transferred mass needs 
to be accreted by the compact star. Part of it may leave 
the system before releasing much potential energy, e.g. via 
a wind from the outer regions of the accretion disk. 

From this it is clear that computing a reliable value for 
r\ is a formidable if not nearly impossible task. Therefore, 
we treat rj as a free parameter and our goal must be to 
arrive at conclusions which are as independent of r\ as 
possible. 

For deriving the criterion for adiabatic stability against 
irradiation-induced mass transfer we follow Ritter (1988) 
(see also Ritter 1996). The only difference is that we use 
now Eq. (JgJ) instead of (|l|) (the latter corresponding to 
Eq. (12) in Ritter's (1988) paper). Accordingly the crite- 
rion for adiabatic stability becomes 



Cs - Cr > Cirr 



= -Ms 



d 



dlnR, 
dM s V dt 
d_ fd\nR t 
dL 



dt 



dL 
ML 



(25) 



The dimensionless number £; rr which is defined by (25) 
measures how sensitively the stellar radius changes in re- 
sponse to irradiation produced by mass transfer. Ignoring 
the influence of irradiation, i.e. setting dL/dM s = 0, gives 
Cirr = and ( p5j ) reduces to the usual criterion for adia- 
batic stability. 

We compute the derivative (d 2 \nR s /dtdL) in ( p5| ) in 
the framework of the bipolytrope model (e.g. KR92) from 
which we obtain (cf Eq. 20) 



d_ 

dL 



dt 



Rs 



GM1 



(26) 



For calculating dL/dM s we shall make use of what we 
shall refer to as the weak irradiation assumption. Mak- 
ing this assumption is tantamount to assuming that at 
any point ($, if) lateral heat transport is negligible com- 
pared to radial transport. Lateral heat transport occurs in 
the form of radiative diffusion and advection because of 
non-vanishing lateral temperature gradients dT/dd and 
dT/dip, and departures from strict hydrostatic equilib- 
rium. We shall show in the Appendix that the weak ir- 
radiation assumption can be justified in those cases we 
are interested in and that we may safely neglect lateral 
heat transport. Accordingly, energy conservation requires 
that at any point (■#, <p) the stellar flux, i.e. the energy lost 
by the star from its interior per unit time and unit surface 
area is 



(27) 



where Ti rr is the effective temperature of the surface ele- 
ment in question. With ( p7|) the stellar luminosity, i.e. the 
energy loss per unit time from the interior becomes 



L = Ri 



2 IT 



F(d,ip) smM'd d(p 



(28) 



Before working out dL / dM s in ([25| ) , we shall first examine 
the reaction of the stellar surface to irradiation in rather 
general terms. For that it is convenient to introduce the 
dimensionless irradiating flux 



x 



F 

± irr 



and the dimensionless stellar flux 
(r = — ! I - x = G(x) 



F 



T ilY (x) 
To 



(29) 



(30) 



where To = T cr f (Fin- = 0) is the effective temperature and 
Fq = F(Firr = 0) = ctTq the stellar flux in the absence 
of irradiation. We note that G(0) = 1 and that we expect 
G(oo) — 0, i.e. that for very high irradiating fluxes energy 
outflow from the stellar interior is totally blocked. Next 
we introduce the function 



dG 

9{X) = "ta 



dF 
dF~ 



(31) 
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which has the following notable properties: First 

g{x)dx = G(0) - G(oo) = 1 (32) 



(33) 



Second, for positive albedos 
< g(x) < 1 Vx>0 . 



Third, if F is a monotonically decreasing function of x 
(and there is no physical reason why this should not be 
so), then 

g'(x) < Vx > . (34) 

From ([52]) and ( [34] ) we can immediately prove through 
integration of g{x) by parts that 



Max [xg(x)\ < 1 



The reason why xg(x) attains a maximum can be under- 
stood as follows: Rewriting xg{x) in dimensional form (us- 
ing (H) and (|l|)) we see that 



xg{x) 



(36) 



F dFir 

The second factor in ([}6|) describes the incremental block- 
ing of the energy loss from the interior with changing ir- 
radiating flux _F; rr . The maximum of xg{x) arises because 
—dF/dFi vr vanishes for large This, in turn, is a conse- 
quence of the fact that as long as the star keeps a negative 
temperature gradient dT/dr in its subphotospheric layer, 
i.e. the superadiabatic convection zone, irradiation can not 
block more than the total flux Fq. 

As we shall see below, the fact that xg(x) has a maxi- 
mum is very important for the stability discussion. In fact, 
we shall see in the next section 5 that for realistic situa- 
tions the maximum of xg{x) is smaller by about a factor 
of two than the strict upper limit given by (|35f) . 

We can now return to Eq. ( p8|) and compute dL/dM s . 
This can be written as 

r-2-K />7T 



dL 

dM, 



-R. 



dF dF, 



dF iri d{-M s ) 



sin'&d'&dip, 



(37) 



where we note that the first factor in the integral is equal 
to -g. With 

dF irr Ftefrtp) 

(38a) 



d(-M s ) {-M s ) 



Air R c a 2 

from ( pi]) we have (using (38a)) 

2 



(386) 



dL 



1 Lo 



Rj, 



dM s 4tt (_M s ) \Ro, 



(39) 



where Lq — 4itRq<tTq is the luminosity of the star in ther- 
mal equilibrium without irradiation. Combining now Eqs. 



5J), (J26|) and (B9h we can rewrite the stability criterion 



as 



Cs — Cr > Cin 



1 T M „ Tfn , ( Rs 

^(Q,ni) — 

47T T K H V Ro 



27T pTX 



x(d, (p) </(x(i9, if)) sm.'&d'&dip 



A^2(C S -C R )—^- 1 (Q,n 1 )(^ 
t m s \ Ma 



(40) 



> 



1 

2^ 



where 



(35) t Ms 



Ms 

-Ms 



x(&,<p)g(x(&,<p)) smMM(p = I ,(41) 



(42) 



is the mass loss time scale. Although the relations ( |40| ) 
and (|4l"| ) do not show an explicit dependence on rj they 
nevertheless depend on it via x. However, because of the 
fact that xg(x) has a maximum, the integral on the right- 
hand side of (^) and ( f4l|) must have a maximum that is 
smaller than Max(xg(x)). Hence we can state that systems 
which fulfill the condition 



1 



A > — Max 
2tt 



. (43) 



x(i?, if)g{x{'d, iff) sin "ddDdtp 

Jo 

which is independent of r), are definitely stable against 
irradiation-induced runaway mass transfer. 

The reason for rewriting the stability criterion ( [io]) in 
the form of ( ]4l| ) or (|4^) is that in the latter conditions A 
does not depend on irradiation but only on the internal 
structure of the donor star (via Cs, tkh, F(Q, n\), Ro) and 
on the secular evolution model (via tm b > Cr an d Rs)- On 
the other hand, all the information about the irradiation 
model is contained in the expression on the right-hand 
side. 

For discussing the stability of mass transfer in the limit 
of very small irradiating fluxes, i.e. x — > 0, we must use 
(38b) instead of (38a) in ((39|). This results in the following 
stability criterion: 

Cs - Cr > Cin- = TZTT1T — HQ, ni) x 



2tt />7r 



JO 



g{x{$,tpf) h($,tp) siwddfldip 



(44) 



> 7T / / g(x(0,<p))h{&,(p) sinddddip . (45) 
27r Jo Jo 

Because of ( |34| ) we arrive at a necessary and sufficient 
condition for stability 

T>—g(Q)[ [ h(<&,<p) smDdDdip . (46) 
27r Jo Jo 
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As above, we have separated in the conditions ([l5]) and 
( f46| ) factors which do not depend on irradiation (collected 
in r) from those which do (on the right-hand side). Like 
A, r depends only on the internal structure of the donor 
star and the secular evolution model. In essence, Eqs. ( (45| ) 
and (^6j) are conditions on the value of r\ in the sense that 
for a given irradiation model, i.e. given h($,if) and g(x), 
rj must not exceed a certain value if mass transfer is to be 
stable. 

Furthermore, Eqs. (Kg) and (46) show in particular 



that because g(x) is maximal for x = (sec Eq. 34), i.e. 
for i<j rr = 0, systems with an unirradiated donor star are 
the most susceptible to irradiation. In other words, if a 
system is stable at the turn-on of mass transfer it will 
remain so later, unless secular effects diminish the value 

of r. 



4-2. Specific irradiation models 

In the following we describe two specific irradiation mod- 
els, a rather simple one, hereafter referred to as the con- 
stant flux or constant temperature model, and a more re- 
alistic one, hereafter referred to as the point source model. 

4.2.1. The constant flux model 

In this model we assume a fraction s(~ 0.5) of the stellar 
surface to be irradiated by a constant average normal flux 



(F in ) 



rj GM(-M S 
8^ R r a 2 



(47) 



Note that the parameter s, introduced above, and s c ff 
which we have introduced in Sect. 3 are not the same 
quantity. However, s and s c ff are related and the corre- 
sponding relation will be given below. 

Comparison of j47| ) with ( pi]) shows that we may write 
ip) as follows: 



M0> <P) 



1 

2 - 
, 



< < n 



< if < 2tt 
< if < 2vr 



(48) 



where now § is the colatitude with respect to the substel- 
lar point and if the azimuth around the axis joining the 
two stars. If the star is assumed to be spherical, the co- 
latitude i9 max of the "terminator" of the irradiated part of 
the surface and s arc related via 



s = - (1 - costf max ) 



With (|48p and (^ 
conditions (^lj) and 

A > 2s(x)g{(x)) 

and 

T>T]sg({x)) . 



(49) 

and (x) — (F ilr )/F , the stability 
45| ) become respectively 

(50) 



(51) 



secondary 
M. 




Fig. 5. Sketch of the geometry of the point source model 
involving a spherical secondary irradiated by a point 
source at a distance a. Note that Fi rr (i?) is the irradiating 
flux normal to the stellar surface. 

These are the results presented earlier in Ritter, Zhang 
and Kolb (1995, 1996). 

Because of ([l|]) not only is the irradiating normal flux 
constant but also the effective temperature on the irradi- 
ated part (hence the name constant temperature model). 
( fl8l) inserted in (|2§| ) yields the luminosity of the star 



L = 4nR 2 s [s(aT^ - (Fir,)) + (1 - s)aT 4 ] 



(52) 



Comparing this modified Stefan-Boltzmann law with (|9[) 
yields the relation between s and s e g: 



ScS 



crT^ 



(F irr ) 



s[l-G((x))} 



(53) 



4.2.2. The point source model 

In this model, which has already been discussed in some 
detail by KFKR96, we assume the donor star to be irra- 
diated by a point source at the location of the compact 
star. For simplicity we assume the secondary to be spheri- 
cal. The chosen geometry is axisymmetric with respect to 
the axis joining the two stars. Denoting again by d the 
colatitude of a point on the surface of the irradiated star 
with respect to the substellar point (for a sketch of the 
geometry see Fig. 5), h($,if) in (pif) becomes 



where 



cos ■& - f s 



(l-2/ s cos0 + /2) 3 / 2 



,f.s 



R s 



fs 



Mi, 



(54) 



(55) 



is the secondary's radius in units of the orbital separation 
a. Because the secondary fills its critical Roche volume, f s 
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is a function only of the mass ratio M c /M s . The termina- 
tor of the irradiated part of the star is at the colatitude 



log T 



$max = arccos(/ s ) 



(56) 



Inserting ( |54| ) in (|d|) or (flq), the stability criteria become 
respectively 



A > / x[tf) g(x{$)) smtidti = I PS 
Jo 

and 



(57) 



T>n g(x('d))h('&) sintidd . (58) 

Jo 

With (^4|) inserted in (^) , the luminosity of the star is 
L = AnR 2 s aT^ x 

U(l + fs) + \J q ™G{x(§)) sin tfdtfj . (59) 

Comparing (|5^) with we find that the effective frac- 
tion of the stellar surface over which the energy outflow is 
blocked is 



Seff = ^ \ 1 _ /* _ 



(60) 



5. The reaction of the subphotospheric layers to 
irradiation 

From the stability analysis we have carried out in the pre- 
vious section it is clear that we need to know more about 
the functions G(x) or g(x) (cf. Eqs. 30 and 31) if we wish 
to use the stability criteria in a quantitative way. So far 
we know only the properties detailed in Eqs. (j32])- (|35|). 
These, however, are insufficient for our purposes. There- 
fore, in this section we shall first use a simple model to 
derive g(x) explicitly and thereafter discuss results of nu- 
merical calculations. As in section 4 we shall adopt the 
weak irradiation assumption. In this approximation the 
relation between the stellar flux F and the irradiating flux 
-F] rr is a purely local property. 

5.1. A one-zone model for the superadiabatic layer 

For determining F(Fi„) we need now a more detailed 
model of the stellar structure than the one we have as- 
sumed in Sect. 3. Whereas in Sect. 3 we have assumed 
with Kippenhahn and Weigert (1994) that the convec- 
tive envelope remains adiabatic up to the photosphere, 
we shall now relax this assumption and take into account 
the existence of a thin superadiabatic convection zone be- 
low the photosphere where convection itself is ineffective 
as a means of energy transport and energy flows mainly 
via radiative diffusion. Sufficiently deep in the star, where 




<i ,4 



log P 



Fig. 6. Sketch of the temperature pressure stratification 
adopted in the framework of the one-zone model described 
in Sect. 5.1. 



convection is effective, i.e. adiabatic, the thermal and me- 
chanical structure of the envelope remain spherically sym- 
metric to a very good approximation even in the presence 
of anisotropic irradiation at the surface (cf. our discus- 
sion in Sect. 2). It is only the very thin superadiabatic 
layer (with a mass of typically 10" 10 M Q ), where energy 
is mainly transported via radiation, which is strongly af- 
fected by irradiation. It is this property which allows us 
to make the weak irradiation assumption, i.e. to treat the 
effects of irradiation in a local approximation and with the 
following simple model. 

For this we adopt for the moment the temperature- 
pressure stratification shown in Fig. 6. This means that 
we replace the "true" structure by one in which we as- 
sume convection to be adiabatic out to a point where the 
pressure is P = Pb and the temperature T = Tb (where 
the subscript B stands for the base of the superadiabatic 
zone). The superadiabatic convection zone extends from 
the point (Pb,Tb) up to the photosphere where P — P p h 
and T = T e g. Because in this zone convection is ineffective, 
we make the simplifying assumption that energy transport 
is via radiation only, i.e. that V = V ra d, where V ra( j is the 
radiative temperature gradient. Because on an irradiated 
part T e s — Tj rr > To, where To is the effective tempera- 
ture in the absence of irradiation, but Tb is assumed to 
be the same irrespective of irradiation, we see that irra- 
diation, by raising the effective temperature, reduces the 
temperature gradient V and thus the radiative energy loss 
through these layers. The superadiabatic layer works like 
a valve which is open if there is no external irradiation 
and which closes in progression with the irradiating flux. 
As sketched in Fig. 6 we assume for simplicity that P p h 
does not depend on Fi TV , i.e. on T e ff. 
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We can now to derive Ti 1T (F rr ) by making a simple 
one-zone model for the superadiabatic layer. For this we 
shall furthermore assume that the optical depth through 
this layer is large enough for the diffusion approximation 
to hold. Then the radiative flux is given by 



F ra d — F — 



dT 4 



3kq dr 



(61) 



where g is the density, k the opacity and the other symbols 
have their usual meaning. 

Now we consider the transported flux in two different 
superadiabatic layers (using subscripts 1 and 2) on an ar- 
bitrary isobar at some level P p h < Pi = Pi < Pb- The 
radiative fluxes are then 



P 



ac ( dT 



3KiQi V dr 
Therefore 

KiQl = F2 (dT i /dr) 1 
n 2 g 2 Fi (dT*/dr) 2 



1,2 



(62) 



(63) 



Now we use for the opacity a power-law approximation of 
the form (nm and the ideal gas equation 



(64) 



P^gT 
I" 



where 1Z is the gas constant and [i the mean molecu- 
lar weight. Since in the regions of interest the dominant 
species (H and He) are neutral, we may assume [i\ — \i 2 . 
Thus with Pi = P 2 we have 

QiTi = g 2 T 2 . (65) 
Inserting now (fLOl) and (|65|) into (|63]) we obtain 



F\ _ T±- b {dT i /dr) l 
~P 2 ~ Tl- h {dT^ldr) 2 



(dT-"/dr) 1 K_ h Ur) 

(dlnT/dr)! r l n 

(d In T/dr) 2 ' — u 



(66) 



Now we make the one-zone approximation by writing 

dT" _ T n (P=P B )~T n (P=P ph ) _ Tg-T c " ff 



dr 



At 



At 



(67) 



d\nT _ In T(P=P B )-\n T(P=P ph ) _ InT B -lnT off 
dr At At 



If we now identify layer 1 with the unirradiated one, i.e. 
set Fl = F and Toff.i = T , and layer 2 with an irradiated 
one, i.e. set F 2 = F and T c g 2 = T'm we obtain by inserting 
@ into (H) 



Fn 1 



F G dFt,. - F 



aT 4 



In T B -In To l _ r. = n 
In T B -In T irr ' " U U 



,(68) 



where we have assumed for simplicity Ari — Ar 2 . Equa- 
tion ( |u8l ) can be solved for T; rr = T ilT (T , T B , F irr ), thus 



providing G(x,Tb). Thus ( |6^) t ogether with ( |28| ) or spe- 
cial cases thereof (Eqs. |2| or \i>9\) can be used as an outer 
boundary condition for numerical computations. 

Let us now briefly consider the consequences of our 
above assumptions that P P h,i = P p ii,2 and Ari = Ar 2 . 
Because the dominant opacity source is H~ and there- 
fore the opacity increases steeply with temperature, the 
photospheric pressure on an irradiated, hotter part of 
the star will be lower than on an unirradiated part (be- 
cause F p h ~ g/ K, where g is the surface gravity). Further- 
more, the hotter the surface, the further out the photo- 
spheric point will be, i.e. Ar 2 > Ari m the above cal- 
culation. Because P p h,2 < ^ph,i> the temperature gra- 
dient V2 below an irradiated part will be higher when 
setting Fph, 2 = Pph,i rather than using the proper value 
of F p h, 2- Therefore by assuming F p h,2 = Pph,i we over- 
estimate the radiative flux on the irradiated part. We 
obtain the same result from using Ar2 = Ari: because 
in reality Ar 2 > Ari, the temperature gradient ( |67| ) on 
an irradiated part is lower than what our estimate with 
Ar2 = Ari yields. Therefore, our very simple one-zone 
model, i.e. Eq. ( |68| ) underestimates the blocking effect 
somewhat. Since, on the one hand, this deficit can eas- 
ily be compensated for by slightly increasing the value of 
b, and since, on the other hand, the precise value of b 
which is appropriate is not exactly determined within our 
model (we shall later take an average value determined 
from published opacity tables), we consider (|6|) a fair ap- 
proximation of the physical situation described. The really 
important aspect of our model is, however, that qualita- 
tively it yields the correct behaviour of a stellar surface 
exposed to external irradiation, and that it is still simple 
enough to allow insight in the situation described. 

We can also compute g = —dG/dx. Differentiation of 
yields 



T o+ 47 r„ ln ( T B/ T °) 



5-6^0 



5-6 = 



(69) 



if T irr < Tb- For consistency with Eq. ( |32j ) we require 
g = if T irr > Tb. 

As can be seen from Eqs. (^8|) and ( |69| ) the functions 
G and g depend only on two parameters characterizing 
the unirradiated star, namely on T and Tb, and on the 
opacity law via b. While To is a well-defined quantity, 
Tb is not because in real stars the run of temperature 
T with pressure F is not as simple as the one assumed in 
our simple model (and sketched in Fig. 6). In particular, 
the transition from convective to radiative energy trans- 
port is smooth and does not occur at one particular point 
as we have assumed in our model. Since Tb stands for 
the temperature at which this transition occurs, we de- 
termine Tb by requiring that in a full stellar model the 
ratio F conv /F ra d of convective flux F conv to radiative flux 
F ra( j reaches a prescribed value, say F conv /F ra d = k. This 
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Fig. 7. Temperature Tb of the point in the superadiabatic 
convection zone of a main sequence star F conv /F la d = k, 
or V rad = (k+ 1)V for k = 0.5 (dashed line), k = 1 (full 
line), and k — 2 (dotted line) as a function of stellar mass. 
Tb is measured in units of To, the effective temperature 
of the unirradiated star. 

is equivalent to the condition V ra d = (k + 1)V. Of course 
the choice of k is somewhat arbitrary but a value k ~ 1 
seems a natural choice. Our simple model is an acceptable 
description of the real situation only if for a given model 
Te(fc) is sufficiently insensitive to k. In order to demon- 
strate that this is indeed the case we plot in Fig. 7 the 
run of Tb as a function of the stellar mass of zero-age 
main-sequence stars with Pop. I chemical composition (X 
= 0.70, Z = 0.02) for three different values of k, i.e. k = 1 
(full line), k — 0.5 (dashed line) and k = 2 (dotted line). 
Fig. 7 shows two important properties of low-mass stars: 
The first one is that a significant superadiabatic convec- 
tion zone exists only in stars with a mass M <; 0.65M Q . 
Below M « 0.6Mq the stratification is essentially adia- 
batic up to the photosphere. This means that application 
of our simple one-zone model is restricted to stars in the 
mass range 0.65M Q M ^ 1M Q . The second property 
shown in Fig. 7 is that the run of Tb(M) for different k is 
qualitatively the same for all three values of k. This means 
that as long as a star has a significant superadiabatic con- 
vection zone the transition from convective to radiative 
energy transport occurs in a rather narrow temperature 
interval, thus justifying our simple approach. 

5.2 Stars with a mass M ^ O.6M 

As Fig. 7 shows, Tb ~ To i n stars with a mass 
M 5s 0.6M Q . This means on the one hand that the to- 



tal optical depth between the photosphere (at T = Tq) 
and the point where T = Tb is small, in fact too small 
for our simple one-zone model, which assumes the dif- 
fusion approximation (Eq. |6l| ), to be applicable. On the 
other hand, this means also that in these stars even in 
the photosphere a non-negligible fraction of the flux is 
transported by convection. Therefore, we must ask how 
irradiation changes the transported flux if the top of the 
(adiabatic) convection zone is at low optical depth. Be- 
cause the convective flux is F conv ~ (V — V a( j) 3 ^ 2 and 
the value of V is directly influenced by irradiation, this 
situation cannot be described by a simple model. Rather 
one ought to determine V by solving the full set of equa- 
tions describing convective energy transport, i.e. in the 
simplest case the equations of mixing length theory. Con- 
sidering the uncertainties inherent in current convection 
theories it is not obvious whether it is possible to make 
general statements about the functions g or G. After all it 
is at least conceivable that already a small irradiating flux 
could result in a significant reduction of (V — V a d) (which 
itself is a rather small number because convection is not 
far from adiabatic) and that therefore dF/dF m could at- 
tain a large negative value. However, as the following ar- 
gument shows, there is a limit to how fast F can drop in 
response to increasing F irr . If F drops too strongly this 
results in an effective temperature T m < Tq. This in turn 
means that the temperature gradient must be steeper than 
in the unirradiated star and, therefore, that F > Fg, in 
contradiction to the starting assumption dF/dF lrr < 0. In 
order to avoid this contradiction T m must not decrease 
with increasing Fi rr , i.e. 9Ti rr /9Fi rr > 0, from which we 
immediately recover (^3|), i.e. g < 1. Because g' < (cf. 
34), g is maximal in the limit F; rr — ■> 0. Therefore we need 
to determine g(0) for the stars in question. For this we 
return to the one zone model (Sect. 5.1). From Eq. (69) 
we find that g(0) increases as Tb/Tq decreases. In fact, in 
the limit T B = T , this model yields 5(0) = 1. Thus the 
closer the adiabatic convection zone reaches to the surface 
the more sensitive the star is to irradiation, i.e. the larger 
g. On the other hand, we know that g < 1 in all cases. It 
is therefore plausible that for stars which are almost adi- 
abatic up to the photosphere, i.e. M ^ O.6M , g(Q) ^ 1. 
As we shall see below this is confirmed by numerical com- 
putations. 



5.3. Results of numerical computations 

Numerical computations of g(x) and G{x) have been car- 
ried out by HR97 for low-mass main sequence stars of 
Pop. I chemical composition using a full ID stellar struc- 
ture code (Hameury 1991), where the outer boundary con- 
dition was changed according to Eqs. ( p7| ) and (|23) with 
Firr = const. The results relevant for this paper are shown 
in Fig. 8 and can be summarized as follows: 
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Fig. 8. The values of g(x=0) (see Eq. [51}) for main se- 
quence stars as a function of stellar mass. Results from 
full stellar models from Hameury and Ritter (1997) are 
shown as a full line, results obtained from the one-zone 
model (Eq. |6^) as short-dashed lines. For comparison we 
show also the run of Max[x<?(a;)] (dash-dotted line) and 
Max(ips) (long dashed line) derived from the results of 
Hameury and Ritter (1997). 

a) With Tb(M) as shown in Fig. 7 and an appropriately 
chosen value for b, i.e. b rj 4 — 5, typical for H~ opac- 
ity, the predictions of our simple one-zone model are 
in good agreement with numerical results as long as 
Ti 1T < Tb- In particular, the run of g(0) shown in Fig. 8 
as a full line for the numerical computations and as a 
short-dashed line for the one-zone model show that 
within its range of validity the latter reproduces the 
numerical results quite well. 

b) g(0) w 1 for stars with a mass M 0.5M Q , as we have 
argued above (Sect. 5.2). 

c) Max(xg(x)) ^ 0.5 over the whole mass range of in- 
terest, i.e. 0.1M Q <; M <, 1M Q . This is shown as a 
dash-dotted line in Fig. 8. 

d) In the point source model (cf Sect. 4.2.2) the integral 
Ips in ( |57j ) has a maximum which is smaller than (1 — 
fs) Max(xg(x)). Max(Ips) is shown in Fig. 8 as a long 
dashed line. As can be seen, Max(ips) ~ 0.20. 

6. Secular evolution with irradiation 

We resume the stability discussion of Sect. 4, but now 
with the knowledge of the function g{x) which we have 
gained in Sect. 5. We shall restrict most of the following to 
the constant flux model (Sect. 4.2.1) and the point source 



model (Sect. 4.2.2). In addition, we shall assume that all 
the properties of a binary along a secular evolution with- 
out irradiation, in particular the functions Y and A defined 
respectively in Eqs. (^l|) and (|45|), are known. Among the 
compact binaries we specifically discuss cataclysmic vari- 
ables, and low-mass X-ray binaries. In particular, we wish 
to examine the following questions: a) which systems are 
stable (unstable) at 1) turn-on of mass transfer, or 2) dur- 
ing the secular evolution, against irradiation-induced run- 
away mass transfer, and b) what kind of evolution do sys- 
tems undergo which are unstable? Because part of these 
questions have been dealt with extensively by King (1995), 
KFKR95, KFKR96, KFKR97 and MF98, we shall mainly 
be concerned with those aspects which have not already 
been treated in detail in the above papers. 

6.1. Cataclysmic variables 

In the following we set M c = A/wd and R c = i?wD , where 
Mwd and i?wD are respectively the mass and the radius 
of the accreting white dwarf. Furthermore we shall restrict 
our discussion to CVs where the secondary is a low-mass 
main sequence star. Because this is the case for the vast 
majority of CVs this is not a very strong restriction. For 
determining the values of A and Y we adopt the standard 
evolutionary scheme for CVs, i.e. the model of disrupted 
magnetic braking (e.g. King 1988 for a review) and use 
results of corresponding model calculations by KR92. 

6.1.1. Stability of mass transfer at turn-on 

Because g(x) is maximal for x = 0, systems are most sus- 
ceptible to irradiation at turn-on of mass transfer. The 
relevant stability criterion for the constant flux model is 
( jsil) and for the point source model Eq. (58). As we have 
already pointed out in Sect. 4 these conditions are in fact 
conditions on the efficiency factor rj because everything 
else is basically fixed. Given the masses of both binary 
components, i.e. Mwd and M s , i?wD follows from the 
mass radius relation of white dwarfs, the secondary's ra- 
dius R s from the mass radius relation of main sequence 
stars, or from the evolutionary history. The orbital sepa- 
ration a follows from Roche geometry via R s and the mass 
ratio, Cs and T from the secondary's internal structure, 
and finally fii again from the mass ratio. Thus Y (Eq. 
[i"5| ) is uniquely determined by Mwd and M s . Given the 
function g(x), and in particular g(0), condition ( |5l| ) is one 
for rjs, and condition (Q) one for r\ only. Taking as an 
example g from our one-zone model with b = 4, s = 0.5 
and Tb /To for k = 1 from Fig. 7, we can plot a line for a 
given r] in the M s -Mwd plane along which mass transfer 
is marginally stable. This is shown in Fig. 9 for various 
values of r\. The figure is to be read as follows: a param- 
eter combination (Mwd, M s ) corresponds to a point in 
Fig. 9. If that point lies above the line corresponding to 
a given value of 77, mass transfer is unstable at turn-on. 
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Fig. 9. Contour lines in the M\vt> — M s -plane along which 
Cirr(-Firr = 0) = 1 for different values of 77. Systems above 
(below) a particular line are unstable (stable) against 
irradiation-induced runaway mass transfer at the turn-on 
of mass transfer. 



What we can infer from Fig. 9 is that for typical WD 
masses Mwd i<i 0.5M q , condition ( pi] ) is violated for sur- 
prisingly small values of 77, i.e. 77 ^ 0.1. This is the case, in 
particular, if M s <; 0.6M Q . If we take instead of ( |5l| ) the 
condition of the point source model ([58|), s is no longer a 
parameter. Taking again the same g as above (one-zone 
model with 6 = 4) the result is qualitatively the same. 
The main difference is that the value of 77 necessary for 
marginal stability needs to be larger by about a factor of 
2. 

The trends seen in the curves of Fig. 9 are easily ex- 
plained: they reflect the run of r(AfvvD, M s ) (cf Eq. 46). 
The smaller T the less stable mass transfer. Because T 
is proportional to -Rwd/Mwd and T^ 1 , and ]?wd/-M\vd 
is a steeply decreasing function of Mwd whereas J 7 ^ 1 ~ 
M ce /M s scales roughly as the relative mass of the convec- 
tive envelope (cf. Fig. 1, dashed line) and thus decreases 
strongly with M s , mass transfer is more likely to be unsta- 
ble (stable) the higher (lower) A%d and M s . The other 
factors entering T, i.e. (Cs~ Cr); M s /R s and (a/R s ) 2 are of 
comparatively minor importance. The fact that the curves 
in Fig. 9 are not monotonic results from the steep increase 
of Tb/T with stellar mass near M = 0.65M© (cf. Fig. 7). 

Because g'{x) < 0, systems in which mass transfer is 
stable at turn on will also be stable against irradiation- 
induced runaway mass transfer for any finite value of F{ TT , 
i.e. M s . The opposite, however, is not true: not all systems 



Fig. 10. The functions A (Eq. |T]) (full line), M&x[xg(x)} 
and Max(Tps) according to the one-zone model (dashed 
lines 1 and 2 respectively), and Max[a; i g(a;)] and Max(Ips) 
derived from the results of Hameury and Ritter (1997) 
(dashed lines 3 and 4 respectively) as a function of the 
secondary's mass along a standard evolution of a CV with 
M WD = 1M Q and M s>i = 1M calculated by KR92. 



which are unstable at turn-on will be so for the secular 
mean mass transfer rate. 

6.1.2. Stability of secular evolution 

The appropriate stability criterion is ( pll"| ) in its most gen- 
eral form, ( |50| ) for the constant flux model and (|5^) for 
the point source model. Again, the left-hand side A of 
these criteria is fully determined by the adopted model of 
secular evolution, whereas the value of the corresponding 
right-hand sides still depends on 77. As we have stressed 
in Sect. 4, we can give a sufficient criterion for stability 
of mass transfer because the function xg(x) has a maxi- 
mum. The corresponding criterion is ( f43|) in its most gen- 
eral form, 



A > 2s Max.(xg(x)) 
for the constant flux model, and 
A > Max(7p S ) 



(70) 



(71) 

for the point source model. Conditions ( ]7C| ) and (|7l]) can 
be read off Fig. 10 where we plot A as a function of sec- 
ondary mass M s (full line) along a standard secular evolu- 
tion of a CV with Mwd = IMq and an initial secondary 
mass M s;i = 1M Q . The evolutionary data are taken from 
KR92. Into the same figure we plot M&x(xg(x)) (= 2s 
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Max(xg(x)) for s = 1/2) and Max(Jps) using results by 
HR97 for g(x) (dashed lines). As can be seen the values 
of Max(xg(x)) and Max(/pg) are almost independent of 
secondary mass, i.e. of the orbital period. For CVs above 
the period gap, i.e. with M s <; 0.25M Q , A increases very 
steeply with decreasing secondary mass. This behaviour 
is mainly due to two factors, first to T~ x ~ M ce /M s 
which increases strongly with decreasing mass, and sec- 
ond to tkh which increases strongly towards lower masses 
mainly because of the mass luminosity relation L ~ M 3 
of low-mass main sequence stars. So, what T essentially 
represents is the thermal inertia of the convective enve- 
lope. As can be seen from Fig. 10, the intersection of A 
with Ma,x(xg(x)) or Max(7 PS ) is near M s = O.7M . Be- 
cause A increases so steeply, small changes in Max(xg(x)) 
or Max(Ips) do not yield a significantly different result. 

When a CV approaches the period gap, i.e. the sec- 
ondary becomes fully convective at M s « 0.25M Q , the sys- 
tem detaches and T — > because tm s ~ * 00 ■ When mass 
transfer resumes A is smaller by typically a factor 10 than 
immediately above the gap, reflecting the reduced angular 
momentum loss below the gap and the fact that at turn-on 
R s = Ro (cf. H|). A starts increasing again with further 
decreasing mass because the Kelvin-Helmholtz time tkh 
becomes very long as the secondary becomes degenerate. 

Because A does not depend strongly on the mass of 
the white dwarf we can generalize the result found from 
Fig. 10: In the framework of standard CV evolution, i.e. 
the model of disrupted magnetic braking, CVs are stable 
against irradiation-induced runaway mass transfer when 
the mass of the secondary star is M s < 0.7M Q . If, on the 
other hand, M s <; 0.7M Q , a system can be unstable but 
need not be so, subject to the value of r\. 

We note, however, that by invoking substantial con- 
sequential angular momentum loss (CAML), i.e. angu- 
lar momentum loss which depends explicitly on the mass 
transfer rate (see King and Kolb 1995 for a discussion of 
CV evolution with CAML), the mass range over which 
systems can be unstable is much larger. This has already 
been pointed out by KFKR96 and confirmed in computa- 
tions by MF98. 

As can be seen in Fig. 10 there is also a slight chance 
that CVs immediately after turn-on below the period 
gap are unstable. According to the constant flux model 
(Eq. 72), some CVs could be unstable, according to the 
point source model, which is more realistic but still very 
optimistic, the stability criterion ( fn] ) is violated only 
marginally if at all. Therefore, our conclusion is that CVs 
below the period gap are very probably stable, at least in 
the framework of standard CV evolution. 

Examining the factors which determine the value of A 
(Eq. 41) we see that within a given evolutionary model 
all factors are determined. Cs, tkh and T depend only on 
the mass of the secondary and its evolutionary history, tm s 
on the adopted rate of angular momentum loss, (s — Cr 
and the evolutionary history (via (dlnR s /dt)th in (§)), 



(r on the mass ratio and R s /Rq on the evolutionary his- 
tory (e.g. Stehle et al. 1996). We see also that making A 
smaller (in order to get the instability for lower secondary 
masses) is possible only by either increasing tm, > i>e. low- 
ering the mean mass transfer rate (M s ), or by increasing 
Cr. The former is practically impossible without upsetting 
the standard evolutionary paradigm for CVs, i.e. the pe- 
riod gap model which requires that above the period gap 
(M 3 £ O.25M , P £ 3 h ) tkh/tm s £ 5 (e.g. Ritter 1984, 
King 1988; KR92; Stehle et al. 1996). Increasing Cr, on 
the other hand, is possible only if a system experiences 
significant CAML. 

6.2. Low-mass X-ray binaries 

At first glance one might suspect that in LMXBs irradia- 
tion of the secondary represents a much larger threat for 
the stability of mass transfer than it does in CVs. How- 
ever, for the following reasons this is very probably not 
the case: first we note the observational fact that very few 
of the LMXBs show X-ray eclipses. This has been inter- 
preted as a consequence of the large vertical scale height 
of the X-ray irradiated accretion disk. This in turn allows 
the secondary to stay permanently in the disk's shadow. 
If this is the case, none or at most a small part of the 
secondary's surface is directly exposed to the accretion 
light source. Second, indirect illumination of significant 
parts of the donor (the high latitude regions or part of the 
back side) is ruled out because this would require a very 
extended scattering corona indeed, with a typical size of 
the scattering sphere (at optical depth r ~ 1) of order or 
larger than the donor star. As a consequence one would 
expect X-ray eclipses to occur much more frequently than 
they are actually observed. Third, heat transport by cur- 
rents from the hot, illuminated to the cool parts in the X- 
ray shadow are probably also negligible. This is because in 
cool stars with a deep convective envelope the superadia- 
batic convection zone isolates the interior from the surface. 
It is itself unable to transport significant amounts of heat 
because of its small heat capacity and the fact that the 
thermal time scale is much shorter than the time scale of 
circulation. 

Heat transport by currents caused by hydrostatic dis- 
equilibrium is, however, of importance in stars with a ra- 
diative envelope. Effects of this are e.g. seen in the X-ray 
binaries HZ Her and V1033 Sco (see e.g. Shahbaz et al. 
(2000), and references therein). 

Because neither indirect llumination nor heat advec- 
tion can contribute significantly to the blocking of the 



stellar flux, the integral I on the right-hand side of (41) 
is much smaller than for a comparable CV, despite the 
fact that in a LMXB there is potentially much more en- 
ergy available for irradiating the secondary. The fact that 
xg(x) has a maximum at x ~ 1 is the reason why even 
large irradiating fluxes do not help. Rather, optimal ir- 
radiation is achieved if as large a fraction as possible of 
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the secondary's surface is irradiated with a flux such that 
xg(x) is near its maximum, i.e. if x rs 1. This is clearly 
not the case in LMXBs. Not only is the surface fraction 
which is directly irradiated small, worse, where the surface 
of the secondary is directly exposed to the X-ray source, 
the associated flux is large, i.e. x » 1 unless n is assumed 
to be very small. The latter is very unlikely considering 
the fact that most of the accretion luminosity emerges in 
form of rather hard X-rays. Thus, unless the secular evo- 
lution of LMXBs is totally unlike that of CVs as far as 
the nature of the secondary and the typical mass transfer 
rates are concerned, the value of A is virtually the same 
as for CVs but, as explained above, the right-hand side of 
(|d[) is much smaller than in CVs. Therefore we conclude 
that very probably LMXBs are stable against this type of 
irradiation-induced runaway mass transfer. As has been 
noted by KFKR97, systems in which accretion is inter- 
mittent rather than continuous, i.e. transient LMXBs, are 
even more stable. 

6.3. Evolution of unstable systems 

Let us now discuss briefly the evolution of systems 
CVs) in which the stationary mass transfer given by 
unstable, i.e. systems for which the stability criterion (4l[) 
or a special form thereof (Eqs. ^o| or ^) is violated. From 
the fact that g'(x) < we know that mass transfer is then 
already unstable at turn-on. Therefore, when mass trans- 
fer turns on, the mass transfer rate increases, and because 
( fi"i"| ) is violated, it does not settle at the secular mean 
(—M s ) given by (|^). However, because the thermal relax- 
ation caused by irradiation saturates both in amplitude 
and with time (se e our discussion in Sect. 3, in particular 
Eqs. (20) and (p3|), and Fig. 4) the mass transfer rate does 
not run away without bound. Rather, there is an upper 
limit: with (20) we obtain for the maximum mass transfer 
rate 



Max(-M s ) « (-M a )+ Ms Maxj^jr 

(Cs - Cr) I tkh 

(,s - Qr tkh 



(72) 



Because mass transfer can be unstable against irradiation 
only if M s il 0.7M Q , i.e. when the secondary has a rel- 
atively thin convective envelope and thus T is large (cf. 
Fig. 1), the maximum contribution of irradiation to mass 
transfer can be several times (sT times) the thermal time 
scale mass transfer rate M s /tkh- Thus, for such systems 
Max(-M s ) > (-M a ). 

After having reached the peak mass transfer rate, mass 
transfer cannot continue at that rate. Rather it must de- 
crease with time for two reasons: first, thermal relaxation 
saturates on the time scale given by (p3|), i.e. (dbxR/dt)^ 
decreases on that time scale. Second, because mass trans- 
fer occurs at a rate above the secular mean, the binary sys- 
tem is driven apart, i.e. (d\n Rr/ dt) > (din R s /dt). Both 



effects result eventually in the termination of mass trans- 
fer. The system becomes slightly detached, the secondary, 
in the absence of irradiation, shrinks. However, because of 
the absence of mass transfer the contraction of the system 
due to angular momentum loss is fast enough to catch up, 
so that mass transfer resumes and the cycle repeats again. 
In other words: if a system is unstable at the secular mean 
mass transfer rate it must undergo a limit cycle in which 
phases of enhanced, irradiation-driven mass transfer al- 
ternate with phases of very low or no mass transfer. The 
conditions for the occurrence of mass transfer cycles in 
semi-detached binaries have been investigated in more de- 
tail and using more general principles (in the framework of 
non-linear dynamics) by King (1995), KFKR95, KFKR96 
and KFKR97. Their main result is that mass transfer cy- 
cles driven by radius variation of the secondary can only 
occur if (d\nR s /dt)th + (d\nR s /dt) nuc in (Q) depends ex- 
plicitly on the instantaneous mass transfer rate. The only 
plausible mechanism providing such a dependence is irra- 
diation of the donor star by radiation generated through 
accretion, i.e. the situation we are studying in this paper. 
The necessary criterion for the occurrence of mass transfer 
cycles found by these authors is identical to what we have 
found here, namely the violation of (pLt|). In the frame- 
work of a linear stability analysis of mass transfer, with 
which we, King (1995), KFKR95, KFKR96 and KFKR97 
have been concerned so far, we can not calculate the long- 
term evolution over time scales tm 3 of systems under the 
irradiation instability. For this the full set of equations de- 
scribing mass transfer and stellar structure under irradi- 
ation have to be solved. Results of such calculations and 
computational details will be presented in the following 
section. Because similar computations have already been 
done by HR97 and MF98, we shall concentrate here on 
aspects which have not been dealt with in detail by HR97 
and MF98, but are, in our opinion, important for better 
understanding of the evolution under the irradiation in- 
stability. 

7. Secular evolution with irradiation: numerical 
computations 

7.1. Computational techniques 

To compute the secular evolution of a compact binary 
with a low-mass donor star we have used the bipolytrope 
programme described in detail in KR92. With respect to 
the procedure described in KR92 we have, however, im- 
plemented two modifications in order to allow for a proper 
treatment of irradiation and its consequences for the sec- 
ular evolution. 

First, we compute the mass transfer rate explicitly 
rather than by using Eq.(||) as in KR92, i.e. we adopt 
the following prescription (e.g. Ritter 1988): 



M s = Mq exp 



Rr — R s 



(73) 
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Expressions for the photospheric pressure scale height Hp 
and the factor Mq in terms of stellar parameters of the 
donor star and binary parameters are also given in Rit- 
ter (1988). In the context of this paper the donor star 
is always a low-mass main sequence star. For such stars 
both, Hp and Mo are only weak functions of the stel- 
lar mass and have typical values Hp/R s rj 1CP 4 and 
Mo ~ lO _8 M0yr _1 . The reason for using ([73]) instead 
of (||) is that the latter is valid only for stationary mass 
transfer. However, this is not a good approximation when 
a system evolves through mass transfer cycles, because 
such cycles proceed unavoidably through phases of non- 
stationary mass transfer. In order to simplify the com- 
putation of —M s we used fixed values for Hp and Mq, 
i.e. Hp = 10" 4 i? s and Mq = Kr 8 M©yr -1 . Further- 
more we used ([73]) also when — M s > Mq, i.e. when the 
donor star overfills its critical Roche volume. (73) is still 



priate choice of h($,<p), i.e. Eqs. (48) or (54), and <?(•$,</?) 
from (J6E 



a reasonable approximation if R s — Rp < few Hp, i.e. 
-M s <> lO^Mgyr" 1 (see e.g. Kolb and Ritter 1990) 
which is adequate for our purposes. 

When computing the mass transfer rate from ( f73| ) we 
use for Hp and Mq, both of which depend on the effec- 
tive temperature, the value of T e ff on the unirradiated 
part of the star, i.e. T c ff = To. There are two main argu- 
ments for this choice: first, in disk-accreting systems the 
inner Lagrangian point L\ is in the shadow of the disk and 
the cooling time of gas in the superadiabatic layer flowing 
from the irradiated parts towards L\ is short compared 
to the flow time in the shadow region. Second, if we use 
T e ff = T; rr instead of T c g = Tq, Hp and —M s would react 
practically instantaneously and with a large amplitude to 
irradiation giving, in turn, rise to a runaway of the mass 
transfer rate on an extremely short time scale (order of 
hours). This is obviously not what happens in the systems 
we do observe. 

Second we use either Eq. ( |52| ) in the constant flux 
model or Eq. ( p9| ) in the point source model of irradiation 
in place of the usual Stefan-Boltzmann law as one of the 
outer boundary conditions for the donor star model. For 
the calculations presented below we have used the one- 
zone model described in Sect. 5 rather than results ob- 
tained from full stellar models by HR97 described earlier. 
Specifically, T irr in (52) or (H|) was computed by solving 
( |68| ) in which Tj rr is determined for an adopted value of 
T] from ( pi|) with h("d, tp) taken respectively from ( |48"1) and 
( [49] ) in the case of the constant flux model, and from ( ]54"|) 
in the case of the point source model. Tp in ( |68| ) is taken 
from the numerical results shown in Fig. 7. Specifically, 
we have used Tp(M s ) for k = 1. Furthermore, in order to 
determine T lrr we need to specify b, i.e. the temperature 
exponent of the photospheric opacity law ([l0|). For most 
of our experiments we have used b = 4 which is adequate 
for H~ opacity. 

The diagnostic quantities Cirr and s c g are calculated 
from (44) and respectively (|53| ) or (|6C| ) with the appro- 



7.2. Computational limitations 

By using the bipolytrope approximation for describing the 
internal structure of the donor star, we can only deal 
with chemically homogeneous stars which have a radia- 
tive core and a convective envelope or are fully convec- 
tive, i.e. with low-mass (M s ^ 1M Q ) zero age main se- 
quence stars. Thus, we are unable to address chemically 
evolved stars, in particular subgiants. Giants, on the other 
hand, have been dealt with approximately by KFKR97. 
We note also that chemically evolved donors among CVs 
might be more common than hitherto thought (Kolb and 
Baraffe 2000; Ritter 2000; Baraffe and Kolb 2000) and 
therefore deserve further study. For our computations, we 
have adopted Pop. I chemical composition. In the context 
of the bipolytrope approximation, the chemical composi- 
tion is of relevance mainly for determining the appropriate 
gauge functions, i.e. polytropic index n% of the radiative 
core and entropy jump h in the surface layers as a function 
of stellar mass (for details see KR92). Even when adopt- 
ing the appropriate gauge functions, we know that the 
value of £s computed in the bipolytrope approximation is 
smaller than what one would obtain for a full stellar model 
if M 0.6M Q . Because of this, in the bipolytrope ap- 
proximation binary systems are more likely to be unstable 
against irradiation than they are in reality if M s ^ 0.6M Q . 

We have used the one-zone model described in Sect. 5. 
As we have discussed there, application of this model is 
practically restricted to stars with a mass M s 0.6 — 
0.7M Q . We have, however, carried out a few calculations of 
systems with a smaller secondary mass. In those cases we 
have assumed g(x) = lif0<x<l and g = otherwise. 
This corresponds to what (^) yields in the limit Tp — > 
To and approximates results of numerical calculations, at 
least for small fluxes, i.e. x < 1, reasonably well. 

Our calculations contain a number of free or at best 
not well-constrained parameters. To mention just the most 
important ones: in the constant flux model we already 
make a very simplifying assumption about the function 
h(i!}, tp), i.e. the irradiation geometry. This assumption re- 
sults in the free parameter s. In addition, we have the 
efficiency factor r\, a parameter about which we know lit- 
tle beyond the fact that probably < 77 1. Using our 
one-zone model (Sect. 5) introduces furthermore the pa- 
rameters Tp and b, both of which can, however, be fixed 
reasonably well by comparison with full stellar models. In 
the point source model wc do not need to specify s. But 
all the other parameters, i.e. rj, Tp and b remain in the 
problem. By using numerical results for the functions G(x) 
and g(x) would get rid of the parameters Tp and b. But 
we would still be left with specifying h("d, tp) and n. On 
top of all that we have also a number of input parameters 
and functions which are already needed for computing a 
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Table 2. Parameters and model assumptions used for 
calculating the evolution shown in Figs. 11-13. Parame- 
ters common to all three examples are: Mwd = IMq, 
i? W D = 510 8 cm, H p /R = 1(T 4 , M = lO^Mgyr" 1 ; ir- 
radiation model: one-zone model with 6 = 4; irradiation 
geometry: constant flux over surface fraction s = 0.5. 



parameter 


Fig. 11 


value 
Fig. 12 


Fig. 13 


M S)i (M Q ) 
V 

rj(yr) 


0.80 
0.035 
1.27 10 8 


0.70 

0.035 

1.2710 s 


0.25 
0.20 
2.95 10 9 



secular evolution without irradiation. Apart from the bi- 
nary's initial parameters the most important of those are 
the angular momentum loss rate and parameters arising 
from assumptions about mass and consequential angular 
momentum loss from the binary system. 

The purpose of the following computations, going be- 
yond the linear stability analysis, is to illustrate the tem- 
poral evolution under the irradiation instability and a 
number of its specific properties which we have uncovered 
in the foregoing discussion. 

As we have demonstrated in the previous section the 
irradiation instability is more likely of relevance for CVs 
than for LMXBs. Consequently, in the examples below we 
have assumed that the accretor is a white dwarf with a 
mass A/wd = 1M and a radius according to the mass 
radius relation (e.g. Nauenberg 1972) of i?wD = 5 10 8 cm. 
We assume that during the secular evolution the mass of 
the white dwarf remains constant, i.e. (Mwd) = 0, and 
that on average the transferred matter leaves the system 
with the specific orbital angular momentum of the white 
dwarf. Other parameters characterizing the three exam- 
ples which we shall discuss subsequently in detail are listed 
in Table 2. 

7.3. Numerical results 

For the first example of a secular evolution with irra- 
diation, results of which are shown in Fig. 11, we have 
adopted the constant flux model with s = 0.5, an initial 
secondary mass M s ^ = 0.8M Q and loss of orbital angu- 
lar momentum on a constant time scale tj = —J/ J = 
1.27 10 8 yr derived from the Verbunt and Zwaan (1981) 
prescription for magnetic braking with / vz = 1 . The value 
rj = 0.035 was chosen such that on the one hand initially 
Cirr > Cs — Cr an d on the other 21 rr < Tp at all times. 
The top panel of Fig. 11 shows the run of the mass trans- 
fer rate with irradiation (full line) and without (dashed 
line), the second the evolution of the secondary's radius 
with irradiation (full line) and without (dashed line). In 
the third panel we show the run of Cirr (full line) and of 
Cs — Cr (dash-dotted line) with time. Finally the bottom 
panel shows the run of s e ff with time. As was to be ex- 
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Fig. 11. Evolution of a CV with an initial secondary mass 
M Si i — 0.8M Q through mass transfer cycles computed by 
adopting the constant flux model with s = 0.5 and the 
one-zone model with b = 4. Further parameters and model 
assumptions are listed in Table 2. Shown as a function of 
time are: in the top frame the mass transfer rate with 
and without taking into account the effect of irradiation 
(full and dashed line respectively); in the second frame the 
secondary's radius with and without taking into account 
the effect of irradiation (full and dashed line respectively); 
in the third frame Cirr according to Eq. ( ffO| ) (full line) and 
Cs — Cr (dashed line); in the bottom frame s e ff according 
to Eq. 

pected for a system in which (initially) Cirr > Cs ~ Cr 
mass transfer evolves through cycles with at least ini- 
tially large amplitude. The amplitude of the radius varia- 
tions, on the other hand, are small, a consequence of the 
small value of Hp/R = 10~ 4 . According to ( |73| ) we have 
A log R s = Hp /R A log M s . We see also that during the 
mass transfer peaks a significant fraction of the stellar lu- 
minosity is blocked. Because rj has been chosen such that 
always T„ r < Tp we always have s c s < s. What is im- 
mediately apparent from this calculation is that though 
the system evolves through mass transfer cycles, these 
are damped on a rather short time scale, i.e. a time scale 
much shorter than tj. The reason for this is that Cirr de- 
creases rapidly with time (mass of the donor). Eventually 
Cirr < Cs — Cr an d mass transfer becomes stable. This is 
mainly a consequence of the increasing thermal inertia of 
the convective envelope, i.e. of a decrease of F/tkh with 
decreasing donor mass. The damping of the oscillations 
could only be overcome if at the same time g increases 
sufficiently strongly with decreasing mass. We shall show 
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below that in a restricted mass range this is indeed possi- 
ble. 

Given the results of the above example we now ask to 
what extent they are representative, i.e. whether the qual- 
itative behaviour depends strongly on the adopted model 
parameters or not. With this end in view we have carried 
out numerous experiments, the results of which we shall 
now discuss. 

Working with the one-zone model (Sect. 5) we ask first 
how the above results change with the parameter 6, i.e. 
the adopted photospheric opacity. Inspection of Eq. ( |69| ) 
shows that for given values of Tb and T irr g(0) increases 
with b. This means that at least for small fluxes the donor 
star is more sensitive to irradiation for larger b. The rea- 
son for this is easy to understand: the larger b the more 
pronounced the increase of the optical depth in the supera- 
diabatic layer in response to irradiation, i.e. of the average 
temperature, and thus the more effective the blocking of 
the energy outflow from the adiabatic interior. Thus, in- 
creasing (decreasing) b above (below) the value b = 4 we 
have used in the example shown in Fig. 11 results in more 
(less) pronounced mass transfer cycles. The time scale on 
which the mass transfer oscillations are damped remains, 
however, essentially uneffected by changes of 6, reflecting 
the fact that .F/tkh does not depend on b. 

Next we compare the results obtained with the con- 
stant flux model (Fig. 11) with those obtained with the 
point source model. The results of a run with the latter 
model and parameters identical to those used for produc- 
ing Fig. 11 (except of s which is not a free parameter in 
this model) are qualitatively very similar to those found 
in Fig. 11. The amount of stellar flux blocked during a 
mass transfer peak is, however, systematically smaller (by 
about a factor of two) than in the constant flux model. 
The main reason for this is that the irradiated area on the 
donor sps = 0.5(1 — f s ) = 0.32 is smaller by about a fac- 
tor of 1.6 than what we have assumed in the constant flux 
example, i.e. s = 0.5. Therefore, in order to achieve the 
same effect as in the constant flux model with s = 0.5, r\ 
in the point source model needs to be increased by about 
a factor of two. Otherwise the results obtained from the 
two models are very similar. Because of this and because 
the constant flux model is computationally much less de- 
manding we have performed most of our simulations with 
that model. 

Next we examine briefly the dependence on the initial 
mass of the donor star. From our extensive discussions 
in Sects. 4 and 6 we know already that below a critical 
mass which is between about 0.6 and 0.7M Q , depending 
on the adopted models, systems following a standard CV 
evolution are stable against irradiation. 

Adopting the constant flux model and the one-zone 
model, a system with an initial secondary mass of 0.7M Q 
is still unstable (whereas with the point source model and 
a more realistic g(x) such a system would be stable). Re- 
sults of an evolution with s = 0.5 and r\ — 0.04 and the 



remaining parameters as in Fig. 1 1 (cf. Table 2) are shown 
in Fig. 12. As can be seen this evolution differs in several 
respects from the one shown in Fig. 11. Initially the ampli- 
tude of the mass transfer cycles increases with time. After 
only a few cycles the peak mass transfer rate is so high 
that the irradiation effect saturates, i.e. T m > Tb and 
Seff = s = 1/2. Although the one-zone model used here 
does not apply when T; rr approaches Tb, the main effect 
of saturation can be modelled anyway by setting g = 
if Tit > Tb- The qualitative behaviour obtained in this 
way remains the same as if a more realistic and smooth 
function g(x) was used. After an initial phase of increas- 
ing amplitudes they later start decreasing and die out very 
rapidly after about 4 10 7 yr. This behaviour can be under- 
stood as follows: We have already pointed out above that 
increasing amplitudes of the cycles can only be expected 
if with decreasing mass g{x) increases fast enough. This is 
exactly what happens in the evolution shown in Fig. 12. 
The fast increase of g(x) is the result of the fast decrease of 
T B /Ti when going from M s w 0.7M Q to M s » O.6M (see 
Fig. 7). Below M s as 0.6M Q , g(x) does no longer change 
much with mass. As a result, the mass transfer cycles are 
then damped because of the secular increase of tkh/-^ 7 - 
When the system eventually becomes stable, irradiation 
is still important in blocking the energy outflow. As can 
be seen from the bottom panel of Fig. 12, s c s is of the 
order 0.37 after the system has stabilized. 

Going to even lower initial secondary masses, mass 
transfer cycles cannot occur unless either £s — or 
(M s ) is much smaller than in a standard CV evolution (cf. 
our discussion in Sect. 6.1). For illustrating this we show 
in Fig. 13 the results of a calculation for which we have 
assumed M s — 0.25M Q , i] — 0.20, s = 0.5 and the much 
smaller angular momentum loss rate of gravitational ra- 
diation. Thus these parameters mimic an unstable system 
just below the period gap. As we have seen at the end of 
Sect. 6.1, CVs just below the period gap can be unstable 
if the constant flux model with s — 0.5 is adopted. They 
are, however, stable if the point source model is used, un- 
less is lowered below the standard value by invoking 
CAML (see MF98). If we use in addition to the constant 
flux model also the one-zone model in the limit Tb — > T\ , 
at low fluxes the donor star is even more susceptible to 
irradiation than if a more realistic form of g(x) had been 
used. However, for the purpose of this exercise it does not 
matter whether CVs below the gap are stable or not. This 
example was chosen just to demonstrate that with a low 
enough (— M s ) and depending on the value of 77 a system 
could evolve through mass transfer cycles as we have con- 
cluded from our stability discussion in Sects. 4 and 6. One 
additional property of this run which is worth mention- 
ing is that the cycles are only rather weakly damped. The 
reason for this is that a star with a mass M s < 0.25M Q 
is always fully convective and therefore T = 7/3 remains 
constant. Thus the decrease of £i rr is mainly due to the 
slow increase of tkh with time. 
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Fig. 12. As Fig. 11, but with an initial secondary mass 
M St i — 0.7Mq. Further parameters and model assump- 
tions are listed in Table 2. 

8. Discussion and conclusions 

With the exception of the numerical examples presented 
in the previous section, where we have assumed the donor 
star to be on the ZAMS, we have so far not made any 
explicit assumption about the evolutionary status of the 
secondary. In fact our stability considerations are valid 
for any type of donor star as long as it has a sufficiently 
deep outer convective envelope. In the following we shall 
therefore briefly address the stability of mass transfer and 
the peak mass transfer rate in case of instability if the 
secondary is more or less nuclear evolved. In addition we 
shall briefly discuss the impact of the donor's metallicity 
on the stability of mass transfer. A detailed investigation 
of the irradiation instability in those cases is, however, 
deferred to a subsequent paper. 

8.1. Effects of evolution and metallicity on the stability 
of mass transfer 

Returning to the stability criterion ( |4l| ) we see that the 
left-hand side A depends essentially only on the ratio of 
two time scales, i.e. the thermal time scale of the con- 
vective envelope r co and the mass transfer time scale tm s 
or, respectively, on the time scale on which mass transfer 
is driven ra r = (2/tj + l/r lmc ) _1 . Mass transfer is more 
stable the larger the ratio T ce /T^ T . 

Evolutionary computations such as those shown in Rit- 
ter (1994) show that at a given secondary mass T ce is the 
shorter the more evolved the star. If magnetic braking ac- 
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Fig. 13. As Fig. 11, but with an initial secondary mass 
M Sj i = 0.25Mq and assuming orbital angular momentum 
loss only via gravitational radiation. Further parameters 
and model assumptions are listed in Table 2. 

cording to Verbunt & Zwaan (1981) is invoked for comput- 
ing the angular momentum loss rate, numerical computa- 
tions (e.g. those shown in Ritter 1994, or King, Kolb & 
Burderi 1996) of the mass transfer from a nuclear evolved 
donor show that Td r is the longer the more evolved the 
donor at the onset of mass transfer. Taken together this 
means that mass transfer is the more unstable (less stable) 
the more evolved the donor. As an example, evaluating 
A along the evolutionary tracks shown in Ritter (1994) 
we find that as long as the secondary has not reached 
the terminal age main sequence, mass transfer is stable if 
P & 5 /l — 6 h . In the sequence starting with a 1.2Mq sec- 
ondary at the end of central hydrogen burning, systems 
down to an orbital period P w 2 h may be unstable. 

Main sequence stars of low metallicity are systemat- 
ically hotter, smaller, more luminous and have a thin- 
ner convective envelope than main sequence stars of so- 
lar metallicity. Hence r CG is systematically shorter for low 
metallicity stars. At the same time tm s is the shorter the 
lower the metallicity. Taken together this means that mass 
transfer is systematically more unstable (less stable) the 
lower the metallicity. As an example, using results from 
Stehle et al. (1997), we find that in CVs with a Pop. II 
donor mass transfer could be unstable for orbital periods 
as low as ~ A h . 

The case of systems containing a giant donor and in 
which mass transfer is driven by its nuclear evolution, i.e. 
T dr = T nuc: na s already been dealt with in some detail 
by KFKR97. We shall mention here only that such sys- 
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tems are systematically more unstable against irradiation- 
induced runaway mass transfer than systems with a main 
sequence donor. There are three reasons for this: First, 
for a giant r co is systematically shorter than for a main 
sequence star of the same mass (because of the much larger 
values of R and L), and this despite the fact that M ce /M 
can be quite large. Second, in terms of radius changes due 
to irradiation, giants react much more strongly than main 
sequence stars. (Note that in the framework of the homol- 
ogy model presented in Sect. 3 the effective v for giants 
is -3!) Third, the mass transfer time scale tm b associated 
with nuclear evolution of a giant is always much longer 
than r co over the entire range of core masses of interest, 
i.e. 0.15M Q <> M c & O.5M . For details see KFKR97 and 
Ritter (1999). 

8.2. The peak mass transfer rate 

Our estimate ([72]) of the peak mass transfer rate achieved 
during a mass transfer cycle can be rewritten as 

2s" 

T 



Max(-M s 



-M s ) 



1 



(74) 



In the case of a fully developed instability it is the sec- 
ond term in ( |74| ) which dominates. Therefore, the peak 
mass transfer rate is essentially determined by the rate of 
expansion of the donor, i.e. 



Max(-M s ) 



(75) 



Thus the lower the thermal inertia of the convective enve- 
lope the higher the peak mass transfer rate. As a conse- 
quence, very high peak rates can result if the mass of the 
convective envelope or tkh is small. If this is the case the 
long-term evolution of the corresponding systems could be 
drastically changed. We note e.g. that even if the donor is a 
main sequence star, the peak mass transfer rate can easily 
exceed the value required for maintaining stable nuclear 
burning on the white dwarf, i.e. A/wd <i lO _7 M yr _1 
(e.g. Fujimoto 1982). With a Pop. I ZAMS donor this 
limit is reached if M s > 0.8M Q (see Fig. 11), with a donor 
at the end of central hydrogen burning if P > 10' 1 or 
M s 0.6M©. These values are derived from the evolu- 
tionary calculations discussed in Ritter (1994). 

The consequences of reaching peak mass transfer rates 
equal to or even in excess of the stable nuclear burn- 
ing limit of the white dwarf can be far-reaching. First, 
if steady nuclear burning on the white dwarf is reached, 
such a system will no longer appear as an ordinary CV but 
rather look more like a supersoft X-ray source. Because 
such systems are bright in the EUV only, they become vir- 
tually undetectable in our Galaxy: very few such systems 
have indeed been found (see e.g. Greiner 1996). Second, 
with nuclear burning on the white dwarf, the nuclear lu- 
minosity exceeds the accretion luminosity by a factor of 
~ 10 — 10 2 . With so much more irradiation luminosity 



available effects other than those discussed in this paper 
might also become important, e.g. driving of a strong wind 
from the donor (see van Teeseling and King 1998, King 
and van Teeseling 1998), and which would change the evo- 
lution of such systems altogether. Third, an unavoidable 
consequence of the very high peak mass transfer rates are 
very extended low states during which a system is prac- 
tically detached and thus virtually undetectable. Fourth, 
from the fact that such systems are barely detectable in 
both the high and low state and that the transition time 
between high and low states and vice versa is very short 
(see KFKR97), one is practically forced to the conclu- 
sion that the observed long-period CVs are either stable 
against the irradiation instability or, for some other rea- 
son do not reach peak rates equal to or larger than the 
stable nuclear burning limit. 

8.3. Conclusions 

In this paper we have studied the reaction of low-mass 
stars to anisotropic irradiation and its implications for the 
long-term evolution of compact binaries. For this we have 
shown that the case of anisotropic irradiation in close bi- 
naries is relevant and that spherically symmetric irradia- 
tion is probably not an adequate approximation. We have 
studied the reaction of low-mass stars to anisotropic irra- 
diation by means of simple homology considerations. We 
have shown in the framework of this model that if the 
energy outflow through the surface layers of a low-mass 
main sequence star is blocked over a fraction s e g < 1 
of its surface, it will inflate only modestly, by about a 
factor ~ (1 — s ff)~ al (see Eq. [l7]) in reaching a new 
thermal equilibrium, and that the maximum contribution 
to mass transfer due to thermal relaxation is s c ff times 
the value one obtains for spherically symmetric irradia- 
tion (i.e. s c ff = 1). The duration of the thermal relaxation 
phase is 0.1r co | ln(l — s e ff)| (Eqs. ^ and ^3[), where t cc is 
the thermal time scale of the convective envelope. In ad- 
dition, we have carried out numerical computations of the 
thermal relaxation process of low-mass stars involving full 
stellar models and using the modified Stefan-Boltzmann 
law (^) as an outer boundary condition. The results of 
these computations (shown in Figs. 2-4) fully confirm re- 
sults from homology and show that the effects caused by 
anisotropic irradiation are not only quantitatively but also 
qualitatively different from those caused by spherically 
symmetric irradiation. 

Next we have carried out a detailed stability analy- 
sis. The criterion for stability against irradiation-induced 
runaway mass transfer in its most general form (arbitary 
irradiation geometry) is given in Eqs. @ and @. One of 
the remarkable results of this stability investigation is that 
it is not arbitrarily large irradiation fluxes which destabi- 
lize a system most effectively. Rather the most effective 
irradiating flux is F m w Fq, where Fq is the surface flux 
of the unperturbed star. 
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The reaction of the stellar surface to irradiation is best 
expressed in terms of a function g{x) — —dF/dFi TT , where 
x = Fin/Fo is the normalized flux. General considerations 
show that Max(xg(x)) < 1. For determining g(x) we used 
a simple, analytic one-zone model for the superadiabatic 
convection zone. The results of this simple model (Eq. ^9|) 
are fou pd to bo in good qualitative and satisfactory quan 
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Appendix A: Lateral heat diffusion 



On the other hand, the radial temperature stratification is 
given by the Eddington approximation (e.g. Tout et al. 1989) 



T 4 (r,tf) 



+ 
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FirrW 



FirrW 



where r is the optical depth. Therefore 
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From the definition of the optical depth we have 
dr 



dr 



= —KQ 
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Inserting (A.5) into (A.4) yields the minimum radial tempera- 
ture gradient in the presence of irradiation: 



Min|(VT 4 ) r | = - ( T 4 r (tf = 0) - Firr(,? ~ 0) 
4 \ a 



HQ 



(A.6) 



Combining now (A.l) and (A.6) we obtain 
(VT 4 )tf . 



Max | 



(VT 4 ),- 1 
4 T 4 r (tf = 0) 



3tf n 



T 4 r (tf = 0) - F irr (tf = 0)/cr kqR s 
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Now, because 4/3i? m ax 
kqHp ~ 1, we have 



Max 



(VT 4 ) r 1 



1 and in the subphotospheric layers 
T 4 r (tf = 0) 



T 4 



Hr_ ... 

R s T 4^ = Q) _ F ilT (<>=0) 

Hp G(x(ti = 0)) + x(tf = 0) - 1 
R s G{x{ti = 0)) 



(A.8) 



Because in the stars in question Hp/R s ~ 10~ 4 is so small, 
the lateral temperature gradient is always much smaller than 
the radial one unless x(i? = 0) 2> 1 and thus G{x($ — 0)) 
<C 1. However, in the range of interest of x and G(x), i.e. 
wher e xg(x) is near its maximum and x w 1 to a few, 
(A.8) yields that the lateral temperature gradient is of order 



(VT)„ « H P /R s {VT) r < (VT) r . 



Here we show that in a realistic situation the lateral temper- 
ature gradient caused by irradiating the secondary anisotropi- 
cally is negligible compared to the radial temperature gradient 
in the subphotospheric layers and that, therefore, lateral en- 
ergy transport by radiation is unimportant. We demonstrate 
this in the framework of the point source model. In that model 
we estimate the lateral gradient of 7] 4 r to be 



(VT 4 r )* 



r 4 r (i9 = o)-r 4 r (i9 = i9 max ) 
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Jl 4 r Q? = Q)-r ( 4 

m ax 
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